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Abstract —We present a simple yet effective routing strategy 
inspired by coverage control, which delays the onset of congestion 
on traffic networks, by introducing a control parameter. The 
routing algorithm allows a trade-off between the congestion level 
and the distance to the destination. Numerical verification of 
the strategy is provided on a number of representative examples 
in SUMO, a well known micro agent simulator used for the 
analysis of traffic networks. We find that it is crucial in many 
cases to tune the given control parameters to some optimal value 
in order to reduce congestion in the most effective way. The 
effects of different network structural properties are connected 
to the level of congestion and the optimal range for setting the 
control parameters. 

I. Introduction 

The ability of traffic networks to support an increasing 
amount of vehicular traffic is becoming crucially important, 
as there are social, environmental and economic consequences 
of poorly managed networks [ 11-[3|. In particular, vehicular 
congestion significantly jeopardizes the performance of traffic 
networks. For instance, there is evidence to suggest that in 
cities such as London, travel delays during busy periods 
increase journey times by an average of 34% 0. Also, a recent 
analysis of road speeds in New York suggests that it suffers 
from a rush hour which lasts all day [5]. There are a number of 
factors that may contribute to congestion on a traffic network. 
The most obvious is an increase in the number of vehicles 
using the network exceeding its capacity, for example during 
peak periods (6). If some roads become congested, congestion 
will tend to spread to other parts of the network, so the faster 
those roads can be cleared, the better the network performs. 
With car ownership in the UK having increased over the past 
40 years [7|, this problem is liable to worsen. 

The advent of automated vehicles offers unique opportuni¬ 
ties to reduce congestion, as control of the route is removed 
from human hands. Improving vehicle routes, according to 
some objective function, could become a new and powerful 
control action that can be taken. The dynamic vehicle routing 
problem has been concerned with this for some time GD-0 
but there are few clearly defined decentralised and adaptive 
control strategies for vehicles. 

In a vehicular network the routing problem is to find 
strategies to assign routes to vehicles in order to minimise their 
travel time and reduce congestion on the network. Strategies 
can be local, if they only rely on information available to the 
vehicle in a local neighbourhood of the network, or global if 
vehicles are assumed to know the state of the entire network. 
The simplest and most obvious global strategy is Dijkstra’s 


algorithm 0. whereby each vehicle takes the shortest path 
between origin and destination, and its extensions (see for 
example p2|-|T5|). All of these modern routing algorithms 
enable extremely fast calculation of shortest paths (down to 
milliseconds), but all of them require some pre-processing 
of the traffic network and storage of data that is produced 
from this pre-processing. They are excellent for large static 
networks, but when a network is dynamic, such approaches 
are unpractical for large networks, and require a high commu¬ 
nications overhead when computed by a central controller and 
communicated to every vehicle. 

Developments have been made for truly novel routing 
strategies. For instance in fT6| , Gzt algorithms were de¬ 
veloped for vehicle navigation, which aim to provide highly 
reliable travel time estimations. The algorithm, named ’Hy- 
perstar’, plans routes based on the probability of encountering 
a delay and needing to divert along a different path. In (T8] 
a distributed algorithm is presented which can calculate the 
socially optimum set of traffic flows, given a set of vehicle 
trips, and then route vehicles probabilistically using those 
optimum flows. 

The aim of this paper is to introduce and explore a simpler, 
yet effective, local routing algorithm based on using not only 
shortest path information but also some knowledge of localised 
congestion on different roads. The idea is that as a vehicle 
approaches a junction, it is assigned the next road on its route 
by evaluating a cost function composed of two terms: one 
related to distance from their destination, the other on the 
measured congestion level of roads out of that junction. The 
approach is inspired by coverage control strategies in control 
theory ]T9| and is compared to other routing strategies on a set 
of representative examples. We find that, when compared to 
more traditional strategies, our approach guarantees a lower 
level of congestion and minimal delays. In particular, we 
uncover a subtle relationship between the mean travel time, 
the level of congestion and the weight in the cost function 
determining the relative dominance of the distance term over 
the congestion term. This seems to suggest that optimal tuning 
of the routing control parameters is possible depending on the 
structural properties of the road network and the car density 
on it. 

II. Coverage Based Routing Algorithm 

Coverage control is a method for controlling distributed 
mobile sensor networks, in order that they can deploy mobile 
sensors and provide optimal coverage of a sensory function 
within a bounded region |20| . In order to maximise the flow 
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Fig. 1. (a) The principle behind coverage based routing. A stream of cars want 
to travel between A and B, whilst another group of cars want to travel between 
C and D. Cars travelling between A and B could take one of two possible 
routes to get there (far left image). As a car approaches a possible turning 
point it evaluates the distance of the journey along each potential route (centre 
image), the yellow arrow indicates the longer route, and a respective cost 
function is shown in (b), for some arbitrary distances. The car also evaluates 
the level of congestion on roads it could take (far right image). The red line 
indicates the road of high occupancy, and the green line indicates the road of 
low occupancy. The respective cost function is shown in (c). In this example, 
one route is shorter, but has higher congestion. The car will evaluate the two 
cost functions and combine them to pick the best route. 


in a traffic network, we hypothesize that it would be preferable 
to distribute vehicles as evenly as possible along roads, in order 
to reduce congestion while keeping their travel time minimal. 
This hypothesis is supported by the traffic flow model proposed 
by Lighthill and Whitham [211, which relates traffic density 
(the number of cars per km of road) to traffic flow (the number 
of cars passing a particular point per unit time), as well as by 
further approaches which have used congestion avoidance as a 
factor for calculating optimal routes for vehicles (22), (23) . In 
our routing algorithm, vehicles attempt to balance two sensory 
functions, one which attracts them towards their destination, 
and another which repels them from roads already heavily 
occupied by other vehicles (see Figure [TJ. 

A. Description of Routing Algorithm 

Consider a set of N vehicles, I = {1,2,..., TV}, in a 
bounded region Q C M 2 . The position of the i -th vehicle is 


1 : if car i is approaching junction Vj 

2: Collect local road occupancy data from a nearby intersection controller 

3: for every road ( Vj,Vk ) G Ej 

4: Compute the distance based sensory function using Q 

5: Calculate the occupancy based sensory function using |2j 

6 : Combine to find the overall cost function using |3j 

7: end for 

8 : Take the road ( Vj , v k ) with the lowest cost function by solving 141 

9: end if 


denoted by the vector pi E Q,i E /. The intended final des¬ 
tination of the i-th vehicle is denoted by some point di G Q. 
Vehicles are located on a road network G = {V, E} contained 
in Q, where V = {1, 2, ...,ra} is the set of m vertices which 
represent junctions, and E is the set of edges which represent 
the roads between junctions. The j- th junction is denoted Vj , 
j G V. An arbitrary road is denoted (vj,Vk) G E, where 
Vj and Vk are the parent and child vertices (junctions) of the 
edge (road) respectively. Each road (vj,Vk) has the following 
attributes: 

1) A load [ L(t)] - modelling the number of vehicles 
currently using the road (vj,Vk) 

2) A capacity [C^ ,k ^] - the maximum number of vehi¬ 
cles that can fit onto road (vj,Vk) 

3) An occupancy [r]^ ,k \t) = L ^ j }y ) ] - the percentage 
of space on road (yj,Vk) occupied by vehicles 

When a car, say the i-th vehicle, reaches a junction Vj, it 
can access two types of sensory functions and p^ ,k \ 

for every road (vj,Vk) G Ej, where Ej is a set of all roads 
connected to junction Vj. The first sensory function, </>) , is 

related to the estimated distance that car i will travel if it takes 
a particular road (vj,Vk) G Ej on the way to its destination: 

= 4>(di,Vj,v k ), 0 < <j> ( ?' k) < 1 


The second sensory function, p^ ,k \ is related to the occupancy 
of road ( Vj,Vk ) G Ej at time t : 

p {j ’ k) = p(v {j ’ k \t)), 0<P^<1 


Expressions 0 and 0 are combined via a tuning param¬ 
eter a G [0,1] in order to assign a cost function to each 

road (vj,Vk) available to car i. When car i reaches junction 
Vj, it computes the cost functions for every road (vj,Vk) G Ej 
as, 


jp) = a ^) + (l-«)pW (3) 

Car i at junction Vj will then take the road ( Vj,Vk ) G Ej such 
that 

min (4) 

(Vj,V k )e Ej 

Table [T| summarises the routing algorithm. Note that by tuning 
the control parameter a we can make the vehicle more or less 
sensitive to distance or congestion respectively. 











B. Choice of Sensory Functions 

The sensory function has been chosen as follows 

_ _ V{v k ,di) +V{vj,v k ) _ 

* max^y V{v h di) + max (S)t)GE V(v s ,v t ) 

where V(vk,di) is the network distance (shortest path) be¬ 
tween junction Vk and destination di , and U(vj,vf) is the 
distance of road (vj,Vk)- Additionally, ma xi e yV(vi,di) de¬ 
notes the largest network distance between the arbitrary junc¬ 
tion vi and destination d{. By contrast, max( S)t ) e£ T>(v s , v t ) 
represents the largest distance of the road in E. The basic 
idea behind this definition is that we want to normalize and 
quantify the effect of distance when car i takes road (vj,Vk) 
and then moves to the destination di along the shortest path. To 
allow for the effect of congestion, the sensory function p^^ k ) 
is chosen as 

11 — otherwise 

where ^ is the critical occupancy of road (vj,Vk), at which 
the onset of congestion begins, and cr is a tuning parameter 
determining the sensitivity of when road (■ Vj : Vk ) goes 

beyond its critical occupancy. 

We were motivated to design the cost function in this way 
in order to limit the in-flow of vehicles to intersections with 
long or growing queues, which is observed in traffic flow 
studies such as (24) to reduce the probability of gridlock. A 
graph of the congestion function is shown in Figure 

[Tc| where the chosen values for the function parameters are 
= 0.2 and <r = 10. 

The choice of a is intended to emphasise the strong cost of 
going above , with an exponential decay used so that we 
maintain some distinction between the cost functions for a set 
of roads which have all gone above their critical occupancy. It 
is worth pointing out that the above sensory functions can be 
selected differently and be made dependent on diverse routing 
factors such as travel time, congestion level and road condition. 

III. Numerical Validation 

To validate the strategy we carried out numerical simula¬ 
tions on some representative examples, first in MATLAB (25) 
and then using Simulation of Urban Mobility (SUMO) (26), a 
widely recognised micro agent simulatior for traffic networks. 
We compare the performance of the coverage based routing 
algorithm with others based only on shortest path computation. 
In particular, our algorithm is contrasted with Dijkstra’s and 
a modified Dijkstra algorithm discussed in (27) in the context 
of communication networks that, when more than one shortest 
path is available, takes the one with the next shortest queue of 
data packets. 

A. MATLAB Implementation 

In our MATLAB implementation, each vehicle calculates 
a new position at every time step, based on the average speed 
of the road it is travelling on (average speed is estimated 
using an occupancy-flow relationship from (28)). When a car 
reaches a junction, it recalculates the optimum route using 




(b) 

Lig. 2. Results from modelling in Matlab. (a) Mean travel time plotted against 
car generation rate for coverage based routing when a = 0.4 (red), shortest 
path routing using Dijkstra (cyan) and a modified shortest path algorithm 
(green), (b) Effect of varying a on mean travel time for coverage based routing. 

the routing algorithm being used. Cars are generated in the 
network according to a Poisson process of rate A, we shall 
term the car generation rate. We denote by A the maximal value 
of the car generation rate at which congestion is observed to 
occur (see Figure |2a). 

Each routing algorithm was tested for increasing car gen¬ 
eration rates in order to see the effect on the routing algorithm 
performance. The behaviour of the modified Dijkstra algorithm 
was very similar to coverage based routing when a = 1. Given 
the correct choice of the tuning parameter a , coverage based 
routing outperformed the other routing algorithms (see Figure 
[2] where we show results for a 3x3 grid network, where all 
vehicles make a trip between the bottom-left and top-right 
corners of the grid). To further test this observation in a 
more computationally effective way and on larger networks, 
we move next to verify the algorithm in SUMO. 

B. SUMO Verification 

To investigate the performance of the algorithm and con¬ 
firm the importance of tuning the parameter a , we used 
SUMO. Using the SUMO API (TraCI), several Python mod¬ 
ules were produced in order to allow vehicles to be rerouted 
dynamically at each junction, according to the selected routing 
algorithm. 

For all simulations, journeys were generated randomly 
using the ’randomTrips.py’ script (included within the SUMO 
tools library) and traffic lights were omitted. In these simu¬ 
lations the car generation rate refers to the number of cars 












entering the network at each time-step, which is constant and 
not a poisson process. 

Mean delay was measured as being more informative than 
the travel time in the scenarios of interest. Indeed trips between 
different nodes had different expected travel times, whereas 
delay was directly comparable between all trips. The mean 
delay for a journey between any Origin-Destination (O-D) pair, 
is calculated as the difference between the fastest possible 
travel time for any vehicle between the O-D pair, and the 
actual experienced travel time for a vehicle making the same 
trip during a simulation. 

The key concepts we wished to investigate were: 

• The existence of an optimal region of a , which 
maximises network capacity and reduces congestion. 

• The performance of coverage based routing against 
routing based on the shortest path only. 

• The effect of different road network structures on the 
performance of the algorithm 

1) Network Topologies : We studied the performance of 
coverage based routing on a selection of network topologies, 
in order to better understand how topology affected its perfor¬ 
mance. All roads in the networks are given the same speed 
limit, so that calculation of the shortest path by distance is 
equivalent to calculation of the shortest path by travel time 
(where travel-time is not updated according to real-time traffic 
data) and there are no ’motorway’ shortcuts. The networks are 
undirected, so there are no one-way streets. 

Tests run using SUMO included a 5x5 and 10x10 node 
grid network, a spiderweb network, a random network, and a 
scale-free network (see Figure [3] and Table [IT] for their structural 
properties). 

The two grid networks analysed are inspired by road 
layouts in places such as Manhattan Island in New York 
City (Google Map - http://tinyurl.com/kz32ut8). Any vehicle 
wanting to travel diagonally across the network can choose 
among a large selection of shortest paths, thereby favouring a 
routing algorithm which considers all the possible paths. The 
random network (Figure 3a| was generated by rewiring the 
10x10 grid network, and links were generated using a similar 
approach to the Erdos-Renyi Planar Graph model in (29). 
The random network is unlikely to offer many paths of equal 
distance between destinations, and so vehicles must balance 
the gains from avoiding congestion, with the cost of taking a 
longer path. The spiderweb network (Figure [3b) is a collection 
of ring roads of increasing size, connected by inward ’spokes’. 
This is a simplification of the road structure found in cities 
such as Beijing (Google Map - http://tinyurl.com/ojyexwk). In 
this road layout many vehicles routing using only the shortest 
path will attempt to use the central ring as a shortcut between 
destinations. Finally, the scale-free network is inspired by the 
evidence that a large number of urban road networks are scale- 
free networks m- 

2) An Optimal Range of a : Figures [4] to [8] show the 
result of approximately 500 simulations for each network 
topology. Specifically, each point in the diagram is obtained by 
evaluating the mean delay computed over 3600 time units (1 
hour) while the car generation rate and the control parameter 





(b) Spiderweb network (c) Scale-free network 

Fig. 3. Network topologies used in the SUMO simulations 
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(a) Mean delay against a and car generation rate (blue = 
lower delay, red = higher delay). Dashed lines indicate 
the sections shown in Figure [4b] 



(b) Cut-throughs of Figure [4a] when car generation rate 
is 1.8 (blue), 2.3 (green), and 2.6 (red). 

Fig. 4. 5x5 Grid: Mean delay as a function of a and A 
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(a) Mean delay against a and car generation rate 
(blue = lower delay, red = higher delay). Dashed 
lines indicate the sections shown in Figure [5b] 



(b) Cut-throughs of Figure [4a] when car generation rate 
is 3.4 (blue), 4.4 (green), and 4.8 (red). 

Fig. 5. 10x10 Grid: Mean delay as a function of a and A 


a are fixed at a certain value. We observe that the existence of 
an optimal value for the parameter a is dependent on network 
topology, and also that when it exists it can be robust to the 
scale of the network. 

We consider that the network becomes congested when 
an increase in the car generation rate causes an exponential 
increase in the mean delay (see Table [In] for the optimal 
ranges identified from these plots, associated to the largest 
value for the car generation rate that the network absorbed 
before becoming congested (A)). 


TABLE II. 


Network 

Topology 

Nodes 

Edges 

Mean 

Degree 

Diameter 

5x5 Grid 

25 

40 

3.2 

800 

10x10 Grid 

100 

180 

3.6 

1800 

Random 

100 

180 

3.6 

1364 

Spiderweb 

50 

95 

3.8 

2000 

Scale-free 

48 

58 

2.4 

1371 


TABLE III. Optimal Range of the Parameter a 


Network 

Topology 

Optimal a 

Value(s) 

Max. Car Gen. Rate 

Before Total Congestion 

Shortest Path | Coverage Based 

5x5 Grid 

0.8 to 0.95 

1.5 

2.3 (+50%) 

10x10 Grid 

0.9 

2.0 

4.4 (+120%) 

Random 

0.25 

2.0 

2.8 (+40%) 

Spiderweb 

0.75 

1.6 

2.6 (+60%) 

Scale-free 

< 0.8 

1.6 

Limit not found 



a 

(a) Mean delay against a and car generation rate 
(blue = lower delay, red = higher delay). Dashed 
lines indicate the sections shown in Figure [6b] 


We have capped the mean delay shown at 500 time 
steps, as although we can observe much higher delays in this 
model, we are most interested in characterising the onset of 
congestion. The phase transition to congestion in the network 
is characterised in part (a) of Figures [?] to [8] by the transition 
between the maroon region (high delay, indicating congestion) 
and the dark blue region (low delay, indicating free flow). 

The results indicate that there are indeed optimal choices 
for a, which can significantly impact the performance of the 
coverage based routing algorithm. In all simulations values 
of a = 0 or 1 performed extremely poorly in comparison 
to any other intermediate value, confirming that it is best 
for vehicles to make local choices based on a combination 
of congestion and shortest path arguments. Setting a = 0 
represents routing based entirely on avoiding congestion, and 



(b) Cut-throughs of Figure [4a] when car generation rate 
is 2.2 (blue), 2.6 (green), and 3.2 (red). 

Fig. 6. Random: Mean delay as a function of a and A 
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(a) Mean delay against a and car generation rate (blue 
= lower delay, maroon = higher delay). Dashed lines 
indicate the sections shown in Figure [7b] 



(b) Cut-throughs of Figure [7a] when car generation rate 
is 1.8 (blue), 2.4 (green), and 2.8 (red). 

Fig. 7. Spiderweb: Mean delay as a function of a and A 

so we would expect a poor performance, because the vehicle 
simply avoids congested roads rather than moving towards its 
destination. Setting a = 1 represents routing based entirely on 
the shortest path (and in the case of multiple shortest paths, 
one will be chosen at random). Our results, see Table II, also 
show a clear dependence of the optimal range of the control 
parameter a on the network structure. 

In the 5x5 grid network (Figure [4a]), the optimal range of 
a is clearly visible in the region 0.80 to 0.95. Within this 
range we observe a consistent increase in network capacity 
indicating that in any application of such a cost based routing 
strategy, ad hoc tuning of the control parameter could have 
a huge impact on performance. When the network is scaled 
up to 10x10 nodes, as shown in Figure [5a] we see that the 
optimal value for a is still 0.9, showing that the choice of a is 
robust to the effects of scale in a grid network. In the random 
network, as shown in Figure [6a] we find that a unique value of 
a does not emerge as optimal and several values of a perform 
equally well, and the same was found in the scale-free network 
(Figure [8a]). However, in the spiderweb network (Figure [7a]) an 
optimal value of a emerged at 0.75. This indicates that both 
the existence of a unique value for an optimal a , and its value, 
are dependent upon the structural properties of the network. 

We can conclude that selection of the tuning parameter has 
a strong influence on the performance of the routing algorithm, 
and that this may be a unique value depending on the network 



a 

(a) Mean delay against a and car generation rate (blue = 
lower delay, red = higher delay). Dashed lines indicate 
the sections shown in Figure [8b] 



(b) Cut-throughs of Figure [8a] when car generation rate 
is 1.6 (blue), 2.2 (green), and 3.0 (red). 

Fig. 8. Scale-Free: Mean delay as a function of a and A 


topology. Even in the random network, it was important to 
choose an intermediate value of a rather than setting it to 0 or 
1 for the best performance. This performance difference can be 
further confirmed by taking vertical cut-throughs in Figures [4a] 
to [8a] and comparing the threshold value of the car generation 
rate at which congestion is observed to occur, for different 
routing algorithms and network structures (see Figures [4b] to 
[ 8 §. 

3) Comparison of Routing Algorithms : We have compared 
our routing algorithm against routing using only a shortest path 
calculation, as calculated using Dijkstra’s algorithm. We find 
that coverage based routing outperforms routing using only the 
shortest path, provided the value of a is tuned correctly. 

Figures [9a] to [9d] show shortest path based routing and 
coverage based routing exhibiting similar delays at low car 
generation rates. However, at some intermediate car generation 
rates we see shortest path routing unable to keep the network 
free from congestion, and delays increase by an order of 
magnitude. In contrast, coverage based routing can continue 
to maintain much lower delays up to some other, higher, 
car generation rate. This demonstrates that coverage based 
routing increases the capacity of the network over routing using 
only the shortest path. In terms of the increase in maximum 
car generation rate that coverage based routing exhibits over 
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Fig. 9. Comparison of mean delay when routing using the shortest path (dashed black line) and coverage based routing (solid blue line) when using an optimal 
value of a from Table mu 


shortest path routing, we see around a 50% increase for the 
5x5 grid network, a 120% increase in the 10x10 grid network, 
a 40% increase in the random network, a 60% increase in the 
spiderweb network, and the absence of a definitive maximum 
car generation rate in the scale-free network (at least over the 
range of car generation rates being considered). 

The simulations indicate that significant gains in network 
capacity can be made by using a coverage based approach 
over shortest path routing, however the size of these gains is 
dependent upon network topology. To further confirm this find¬ 
ing and allow a better comparison between different network 
structures, in the final version of this paper we will include 
simulations on a further set of networks such as a scale-free 
rewiring of the 10x10 node grid. 

IV. Conclusion 

The routing strategy we have presented here shows how 
a simple cost function coupling global information about dis¬ 
tances with local road occupancy data can yield improvements 
over shortest path routing. These improvements are an increase 
in the capacity of the network, and hence the ability to avoid 
congestion which leads to delays. We have also shown that 
in order to optimise this routing strategy a control parameter 
must be chosen appropriately, and this value appears to depend 
on the network topology. We envisage that the algorithm could 
be deployed via a junction infrastructure able to communicate 
to each vehicle arriving at the junction the best road to take 
next. We wish to emphasise that the algorithm we present here 
can be deployed effectively and is easily scalable because of 
its simplicity. Future work will address the possibility of each 
car, or junction, becoming able to tune the control parameter 


a in real time via local adaptive strategies aimed at further 
minimising congestion and guarantee fairness. 

Also, the routing strategy will be tested on real road 
network topologies, using realistic (asymmetric) loading on 
the network, including a variety of vehicle types. 
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